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Abstract 

Aii'-Kauffman networks £^ are a subset of the Boolean functions on 
N Boolean variables to themselves, A^v = : Z^-^Z^}. To each NK- 
Kauffman network it is possible to assign a unique Boolean function on 
variables through the function \E' : £^ A^v- The probability Vk that 
\1/ (/) = \E' (/'), when /' is obtained through / by a change of one of its K- 
Boolean functions {bx '■ Z^ —>■ Z2), and/or connections; is calculated. The 
leading term of the asymptotic expansion of Vk, for A ^ 1, turns out to de- 
pend on: the probability to extract the tautology and contradiction Boolean 
functions, and in the average value of the distribution of probability of the 
Boolean functions; the other terms decay as O (1/A). In order to accomplish 
this, a classification of the Boolean functions in terms of what I have called 
their irreducible degree of connectivity is established. The mathematical find- 
ings are discussed in the biological context where, \Ef is used to model the 
genotype-phenotype map. 



Short title: Robustness Against Mutations. 

Keywords: Cellular automata, irreducible connectivity, binary functions, 
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1. Introduction 



TVii'-Kauffman networks are useful models for the study the genotypc- 
phenotype map \E'; which is their main application in this work An NK- 
Kauffman network consists of Boolean variables Siit) e Z2 (i = 1, . . . , A"), 
that evolve deterministically in discretized time t — 0, 1, 2, . . . according to 
Boolean functions on A' (0 < A' < A^) of these variables at the previous 
time t — 1. For every site a A'-Boolean function /,j : Z2 is chosen 

randomly and independently with a bias probability p {0 < p < 1), that 
fi — 1 for each of its possible 2^ arguments. Also, for every site i, K 
inputs (the connections) are randomly selected from a uniform distribution, 
among the A^ Boolean variables of the network, without repetition. Once the 
K inputs and the functions fi have been selected, a Boolean deterministic 
dynamical system; known as a AA-Kauffman network has been constructed. 
The network evolves deterministically, and synchronously in time, according 
to the rules 

S,{t + l) = h {Si, (t), Si,{t), Si,{t)) , i = l,...,N, (1) 

where i,„ 7^ i^^, for all m, n = 1, 2, .... A', with m ^ n; because all the inputs 
are different. An AA-Kauffman network is then a map of the form 

Let us denote by the set of AA-Kauffman networks that might be built 
up, by this procedure, for given A" and K. They constitute a subset of the 
set of all possible Boolean functions on A^-Boolean variables to themselves 

In Ref. 1, a study of the injective properties of the map 

* : £^ ^ Ajv (2) 

was pursuit for the case p — 1/2; where the Boolean functions are extracted 
from a uniform distribution. Using the fact that An = ^2'^^ where is 
the set of functional graphs on 2^ points ^] the average number (A^, K) of 
elements in that ^ maps into the same Boolean function was calculated ^. 
The results showed that for, X ~ O (1) when A" ^ 1, there exists a critical 
average connectivity 

/C=^log.log,(|^)+o(^); (3) 
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such that i^{N,K) e'^^ > 1 (93 > 0) or 'd{N,K) ^ 1, depending on 
whether K < Kc or K > Kc, respectively. 

In genetics, NK-KauSrasm networks are used as models of the genotype- 
phenotype map, represented by \E' The genotype consists in a particular 
wiring, and selection of the Boolean functions fi in (1), which give rise to the 
A'^/f-Kauffman network; while the phenotype is represented by their attrac- 
tors in ^l' (C^) C Ajv — ^2^^ "^"^^ The K connections represent the average 
number of epistatic interactions among the genes, and the Boolean variables 
Si; the expression "1" or inhibition "0" of the i-th gene. 

A well established fact in the theory of natural selection is the so-called 
robustness of the phenotype against mutations in the genotype At 
the level of the genotype, random mutations (by radiation in the environ- 
ment) and recombination by mating, constitute the driving mechanism of 
the Evolution Theory. Experiments in laboratory controlling the amount of 
radiation have shown that; while the change in the phenotype vary from 
species to species, more than 50% of the changes have no effect at all in 
the phenotype In Ref. 1 it was shown that the signature of genetic 

robustness can be seen in the injective properties of with a many-to-one 
map representing a robust phase. This happens if X < with the value of 
N to be substituted on (3), determined by the number of genes that living 
organisms have. This number varies from 6 x 10^ for yeast to less than 4 x 10^ 
in H. sapiens ^. Substitution in (3) gives in both cases that K < 3 ^. 

In this article it is calculated; for a general bias p, the probability Vk that 
two elements /, /' G such that /' is obtained from / by a mutation, 
give rise to the the same phenotype, i.e. ^ (/) = ^ (/'). For a mutation, it 
is intended a change in a Boolean function /j, and/or its connections. The 
results impose restrictions in the values that K and p should have, in order 
that Vk ^ 1/2; in accordance with the experiments. 

The article is organized as follows: In Sec. 2, I set up a mathematical 
formalism that allows to write (1) in a more suitable way for calculations. In 
Sec. 3, the expressions of the different probabilities involved in the calculation 
of Vk are established. In Sec. 4, I introduce a new classification of Boolean 
functions according to its real dependence on their arguments; which I call 
its degree of irreducibility. This classification is used in Sec. 5 to calculate 
the invariance of Boolean functions under changes of their connections and 
so; calculate Vk- In Sec. 6 the conclusions are set up. In the appendix, 
two errata of Ref. 1 are corrected, and it is shown that they do not alter 
the asymptotic results of Ref. 1. So, the biological implications there stated 
remain correct. 
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2. Mathematical Framework 



Now we introduce a mathematical formalism with the scope of write (1) 
in a more suitable notation for counting. All additions between elements of 
1-2 and its cartesian products are modulo 2. 

Let A^iv = {1, 2, . . . , A^} denote the set of the first N natural numbers. A 
X-connection set Ck-, is any subset oi M.n with cardinahty K. Since there 
arc (^) possible A'-connection sets; we count them in some unspecified order, 
and denote them by 

C^^ = {ii, i2,..., iK} C Mn, with a = 1, . . . , (J^^ , (4) 

where, without lost of gcncrahty; ii < 12 < ■ ■ ■ < ix, with 1 < im < A" 
(1 < m < A'). We also denote as 



-pN 
J- K 



k J a=l 



the set of all X-connections. To each X-connection set Cj^^ it is possible to 
associate a X-connection map 



■ "-2 "-2 1 

defined by 

(S) = Ct'' . . . , = . . . , 5,,) V S e 2^. 

Any map 

hK-.TL^ ^ 22, (6) 

defines a A'-Boolean function. Since = 2^; bx is completely determined 
by its X-truth table T^, given by 

Tk = [Ak hK] , 

where is a 2^ x K binary matrix, and is a 2^ dimensional column- 
vector, such that: 

The s-th row (1 < s < 2^) of A^ encodes the binary decomposition of s, 
and represents each one of the possible 2^ arguments of bx in (6). So, Ak 
satisfies ^ 

K 



1=1 
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And 

bif = [(7i, (72, . . . , (72k] , (7) 

where (Jg E '2.2 {1 < s < 2^), represents the images of (6). 

There much as 2^ i^-truth tables corresponding to the to- 

tal possible vectors (7). i^-Boolean functions can be classified according to 
Wolfram's notation by their decimal number jj, given by ^'^^ 



1 + (8) 



Let us add a superscript (/i) to a K-Boolean function, or to its truth table, 

whenever we want to identify them. So, and T^^ refer to the /x-th K- 
Boolean function and its truth table respectively. Within this notation, the 
set of all X-Boolean functions Ek is expressed as 

Of particular importance are the tautology ft^"* = 6^ ^ and contradiction 
^^K = ^^K i^-Boolean functions; with images: 

= [1, 1, . . . , 1] (9a) 

and 

bi,^) = [o,o,...,o]. {%) 

Table 1 gives an example of the A2 matrix (representing the four possible 
entries of 5*1 and 5*2) with the sixteen possible 2-Boolean functions listed 
according to their decimal number (8). 
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Table 1. The A2 matrix, with the sixteen 2-Boolean functions. 



Within the preceding notation, the dynamical rule (1), now may be rewrit- 
ten as r ^ r s 

S, (t+l)^ O Ct'^ (S {t)) , i = 1, . . . , TV; (10) 

where, some of the indexes ccj and //j may be equal for different values of i, 
and S {t) e Z^. 



3. The Invariance of A^if-KaufFman Networks 

Now we are interested in calculate the probability Vk, that (10) remains 

invariant under a change of a connection and/or a X-Boolean function 

h^^ . So, we must study the number of ways in which this could happen; i.e. 
what conditions should prevail in order that for some i, 

h^^Kc*^'''\s) + h^;^Kc'i^'\s) = VSez^, (11a) 

for ctj 7^ /3i and/or Hi ^ Vi. Let us use a shorthand notation and skip to write 
the indexes ctj, and \Xi. Then (11a) may be written as 

Ik o C*k (S) + o (S) = V S e Z^, 

where, Ik = bx + ^bK and = + AC^; with Afe^ e Ek, and AC^ a 
X-connection map. Explicit substitution gives 

bK o AC*K (S) + AbK o C*K (S) ^ VSez^. (11&) 

Equation (lib) could be satisfied in three different ways: 

i) Event A: A change in a X-connection Ck without a change in a K- 
Boolean function bK- This implies AbK = VS e zf", and from (9b) 

^bKoAc*K{s)^b'^;;\ 
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ii) Event B: A change in a A'-Boolcan function without a change in a 
X-connection Ck- This impUes ACjj^- = VS G 2^. From (9b) there 

follows AbK o (S) = VS e ^ AbK = bP. So, the /T-Boolean 
function must remain unchanged. 

Hi) Event C: A change in a K-Boolean function bx and a change in a /T- 
connection Ck- In this case, both ACj^ 7^ and A6;^ 7^ 0; and also 
(lib) must hold with independency of A, and B events. So, from (9a), 

it must happen that bx ° AC^ — b^^p , and Abx o (7^ = 6^"* . 

Since the events A, B, and C, are independent, the probability Vr that 
(11) are satisfied, is given by the combined probabilities P{A), P{B), and 
P (C) that A, B, and C happen. So, 

Vk = P{A)+P{B) + P{C)-P{A) P{B)-P{A) P{C)-P{B) P{C) 
+ P{A) P{B) P{C). (12) 

For a general bias p (0 < p < 1) that = 1, for 1 < s < 2^ in (7), the 
probability IT {bx) to extract the X-Boolean function bx is given by 

n(M=p'^(l-p)'"-" , (13a) 

where 

ci; = a; (6ii:) = ^cTs, (136) 

is the weight of ft^. 

The following considerations are in order: 

i) P {A) is the probability that the projected function 

bt^ = bKoCt^:7L^^7i:^ (14) 

remains invariant under a change of the i^-connection. To get read 
of this we must first introduce the concept of irreducibility of Boolean 
functions, which is going to be done in the next section. 

ii) P (B) is the average probability that bx remains invariant by a muta- 
tion, given that bx has occurred. Then 



Since there are (^^ ) X-Boolean functions with weight a;, from (13) we 
obtain 

P{B)^Y1 D P""" (1 - Pf"'"'' = [1 - 2P (1 - P)f . 

ui=0 V ^ / 

Hi) P (C) is the probabihty of extracting twice the tautology A'"-Boolean 
function. Prom (9a) and (13b) from (13a) 

So, we obtain the asymptotic expression for (12) 

Vk^P {A) + [1 - 2p (1 - p)f [1-P {A)] + O {p'"^') , (15) 
for > 1. 

4. The Irreducibility of the Boolean Functions 

Not all the fC-Boolean functions depend completely on their K arguments. 
For instance, let us consider the 2-Boolean functions of table 1: Rules 1 and 
16 {contradiction and tautology, respectively) do not depend on either 5*1 or 
while rules 6 and 11 {negation and identity, respectively) only depend 
on 5*1. Due to this fact, let us make the following definitions: 

Definition 1 

A X-Boolean function bx is reducible on the m-th argument Sm (1 < < 
K),ii 

bx {Si, . . . , Sm, • • • , Sk) = bx {Si, . . . , Sm + !,•••, Sk) V S G zf . 
Otherwise, the ii'-Boolean function bx is irreducible on the m-th argument 

Sm- 
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Definition 2 

A X-Boolean function bx is irreducible of degree A (0 < A < K); if it is 
irreducible on A arguments and reducible on the remaining K — X arguments. 
If X = K, the iT-Boolean function is irreducible. 

Let us denote by Ik (A) the set of irreducible /^-Boolean functions of 
degree A. From definitions 1 & 2, may be decomposed uniquely in terms 
of X;, (A) by 

K 

Ek^IJ Ik (A) , (16a) 

A=0 



with 



Ik (A) n Ik (A') =0 for X^ A'. 



(166) 



The cardinalities /Sk (A) = ^^Ik (A) may be calculated recursively, noting 
that /9i^(A), must be equal to the number of ways to form A irreducible 
arguments from K arguments. This amounts to (^) times the number of 
irreducible A-Boolean functions /3x (A); thus 



Pk{X) 



/3a (A) 



(17) 



Setting K — X in (16a) and calculating the cardinahties, follows that 



A-1 



2''=E/^^(^)+/^a(A). 



i/=0 



Substituting back into (17) the following recursion formulas for the number 
of irreducible X-Boolean functions of degree A are obtained 



Pk{X) 



A-1 



2^^-E/^a(-) 



1^=0 



(18a) 



and 



{18b) 



Note from (9), that and are irreducible of degree zero. So from 
(18b), 

Ik (0) ^ {bi\ bP} . (19) 
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Some first values for Pk (A) are: 

(3k{2) = 5K{K -1), 
109 

/3k(3) = -^X(X-1) iK-2), 
907 

(4) = ^2-^ - 1) - 2) (i^ - 3) , 

etc. 



5. The Probability P {A) 



Let us now calculate P {A) to obtain Vr from (15). The probability P {A), 

that defined by (14), remains invariant against a change in de- 

pends in the degree of irreducibility of bx', i-e. on which of its K argu- 
ments it really depends. To calculate it, let us first calculate the probability 



P 



A b*^"^ = 0|6x G Tk (A) that, b*^°'' remains invariant due to a change in 
the X-connection Cj^^; given that 6^ is irreducible of degree A. 
Let bx e Ik (A) be irreducible in the arguments with indexes 

nil, 1712, ■ ■ ■ , 'iTT'X , where mi < m2 < • • • < nix 

such that 1 < mi < K {1 < I < X). Let us also rewrite (4) more explicitly 
putting the superscript (a) into its elements; then 

cP = {i^r\i''\---,^P}^M^. 

Now, associated to b*^'^\ we can define its X-irreducible connection by 



Within this notation the set (b^'^^^ , oi the if-connections C]^-* that leave 



b^ ' invariant, is given by 

[bt') = [cf e I zg) = V / = 1, 2, . . . , a} . (20) 
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Then 



Ab 



*{a) 
K 



(a) 



211 



0\bK e Ik (A) 

Prom (5), = (^). To calculate #0^ (&^°^), let us note that the 

i^'-connections C^^ G (^^jc"'*^ have A elements fixed, the elements of 
Jx ^6^"'*^ , and K—\ elements free, which are the elements ofM.N\Jx (^b*^"^^ ■ 
Thus, (b*K^^^ equals the number of subsets of A^iv\ (^b*^°'^^ that can 



be constructed with K — X elements. Since 



MN\Jx{bt^)] =N-X, 



we obtain 



A^- A 
K-\ 



(22) 



That only depends in the degree of irreducibihty A of bx and not in the 
connection index {a). Substituting (22) into (21) we obtain 



A 6 



K 



0|6^ e Jif(A) 
Due to (16), P (A) is given by: 



K\ {N-X)\ 
N\{K-X)\ 



(23a) 



K 

P (^) = J] P [a 6^") = 0\bKe Ik (A)] P [bK e Ik (A)] , 



A=0 



(236) 



where P [bK £ Ik (A)] is the probability that bK be irreducible of degree A. 

The value of P [bK G Ik (A)] depends on Pk (A) [calculated from (18)], as well 
as on the particular way in which the fT-Boolean functions bK are extracted. 

When K O (1) for N ^ 1, equations (23) behave asymptotically like 



p{A)^p[bKeiKm+o\^- 

So from (19), the leading term of P (A) comes from the probability to extract 
the tautology (9a) and contradiction (9b) i^T-Boolean functions. We obtain 
from (13) 

/ 1 



p{A)^p' +{i-py 



o 



N 
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Prom (15) the probability that (1) [or equivalently (10)] remains invariant by 
a change on a i^-Boolean function and/or its connection; is given by 



Vk ~ +{i-pY 

+ [l-2p(l-p)]'"{l- 



+ (1 



-Pf]] 



+ 



(24) 



6. Conclusion 



A classification of X-Boolean functions in terms of its irreducible degree 
of connectivity A was introduced. This allowed us to uniquely decompose 
them through (16), and calculate the asymptotic formula (24) for Vr] that 
an iVif-Kauffman network (1) remains invariant against a change in a K- 
Boolean function and/or its /C-connection. Figure 1 shows the graphs for 
Vk vs p; for different values of the average connectivity K. The graphs 
attain a minimum and are symmetric at p = 1/2 (the case of a uniform 
distribution). For p fixed, Vk < Vk' for K > K' . 

These results are specially important when TVif-Kauffman network are 
used to model the genotype-phenotype map (2) Experiments to study 
the robustness of the genetic material have been done by means of induced 
mutations ^^^i. The results varied among the different organisms studied, 
but it is estimated that in more than 50% of the cases the phenotype appears 
not to be damaged. In A^if-KaufTman networks this phenomena is manifest 
when Vk > 1/2. Figure 1 shows that is possible to be in agreement with 
the experimental data without a bias {p = 1/2), provided K < 1.25 for the 
average connectivity. For the case K = 2 this happens only for values of 
p outside the interval [0.21,0.78]. There is no surprise that biassed values 
of p increment the value of Vk since they tend to increase the amount of 
tautology and contradiction functions (9) through (13). 
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Appendix: Errata in Ref. 1 

All quotations to equations in Ref. 1 are preceded by an "R" , those intro- 
duced here by an "A" , while all the others refer to equations of the present 
article. 

In Ref. 1 it was wrongly stated that the only Boolean functions that con- 
tribute to the number of redundances r in (R16) are: the tautology, the 
contradiction, the identity and the negation. In fact there are contributions 
from many more functions, their number growing with K for K < N [in 
the case K = N oi the random map model r = 0; as explained further); 
according to their classification in terms of its degree of irrcducibility defined 
in Sec. 4 of this article. Furthermore; the contribution to r of the identity 
and negation functions were calculated as 2N [(^ij^) — l] , while the correct 
value is 

Nevertheless these inconveniences: 

• In the asymptotic expansion of (R18) for N ^ 1, the contribution O (1) 
is originated from the tautology and contradiction functions. 

• While the wrong reported contribution 2N — l] , of the identity 
and contradiction fvLiactioiis, turns out to be O (1), it just adds an extra 
term In [Kc + 1) in (R22) that does not contribute to the O (1) term 
of its solution (R23). However it gives a wrong, and slower, decaying 
error O (In In In A^/ In AT). 

• The rest of the Boolean functions, with A > 2, give an O (l/N'^) con- 
tribution to (R18). 

This implies that all the asymptotic results and their genetical conse- 
quences remain correct; while the decaying error term in (R23) becomes 
0{l/N\nN) since the correct value (Al) gives a contribution 0{1/N) to 
(R18). 

The correct results are obtained as follows: 

From (18) and (20), the number of redundances that the elements of (A) 
furnish is given by i3k (A) [#e^ (A) - l] . From (22), the correct value of r 
is: 

K 

r = Y,PK{X) 



(Al) 



A^- A 
K-X 



- 1 



(A2) 
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Note that: 



• The contribution of A = 0, is the one that corresponds to the tautology 
and contradiction X-Boolean functions. 

• The contribution of A = 1, is the one given by (Al), with (1) = 2/^ 
obtained from (18). 

• The contribution of A = is zero. So, irreducible K-Boolean functions 
give raise to injective maps. 

• In the special case of the random map model 3.5,i3. 7- = Q as it should 

be, due to the fact that, for such * : ^ Ajv defined by (R4) 

[respectively by (2) in this article], becomes a bijection so 

where Q2N is the set of functional graphs from 2^ points to themselves ^ . 

With this background, the correct equations (R17), (R18), (R19), (R22), 
(R23), and (R25); are given as follows: 

Prom (R16) and (A2) we obtain 

K 



(z:^) = 2^ 



A=0 



N-X 
K -X 



N 



Now 



{N,K) = {l-^{N,K)f , 
with ip depending also on N; and given by 

Eto/5i.(A)[0-l] 



(R17) 
(R18) 

(R19) 



Prom (18), (p {N, K) admits for A?" ^ 1 the asymptotic expansion 



22^-1 

that gives for the equation 'i?~^(A", Kc) = 1/2, of the critical connectivity, 

1 



2N 
1^ 



1 + 



N 



(R22) 
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The solution of (R22) is now 

^««lofelog.g)+o(^). (R23) 



And (R25) is now given by 



Ai^c^ ., ,,3, ...... (R25) 



1 

(ln2)Mog2(2A^/ln2) ^ \lnN 

This shows that the asymptotic formulas, and conclusions of Ref. 1 are 
correct. 
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Figure caption: 

Figure 1. (Color online), Graphs for Vk vs. p for different values of the 
average connectivity K. = 1 in red, K = 1.25 in green and K — 2 in blue. 
The important Vk = 1/2 value, is in magenta. 
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